home *** CD-ROM | disk | FTP | other *** search
/ EuroCD 3 / EuroCD 3.iso / Programming / Python1.4_Source / Objects / longobject.c < prev    next >
Encoding:
C/C++ Source or Header  |  1998-06-24  |  29.3 KB  |  1,433 lines
  1. /***********************************************************
  2. Copyright 1991-1995 by Stichting Mathematisch Centrum, Amsterdam,
  3. The Netherlands.
  4.  
  5.                         All Rights Reserved
  6.  
  7. Permission to use, copy, modify, and distribute this software and its
  8. documentation for any purpose and without fee is hereby granted,
  9. provided that the above copyright notice appear in all copies and that
  10. both that copyright notice and this permission notice appear in
  11. supporting documentation, and that the names of Stichting Mathematisch
  12. Centrum or CWI or Corporation for National Research Initiatives or
  13. CNRI not be used in advertising or publicity pertaining to
  14. distribution of the software without specific, written prior
  15. permission.
  16.  
  17. While CWI is the initial source for this software, a modified version
  18. is made available by the Corporation for National Research Initiatives
  19. (CNRI) at the Internet address ftp://ftp.python.org.
  20.  
  21. STICHTING MATHEMATISCH CENTRUM AND CNRI DISCLAIM ALL WARRANTIES WITH
  22. REGARD TO THIS SOFTWARE, INCLUDING ALL IMPLIED WARRANTIES OF
  23. MERCHANTABILITY AND FITNESS, IN NO EVENT SHALL STICHTING MATHEMATISCH
  24. CENTRUM OR CNRI BE LIABLE FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL
  25. DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR
  26. PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER
  27. TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR
  28. PERFORMANCE OF THIS SOFTWARE.
  29.  
  30. ******************************************************************/
  31.  
  32. /* Long (arbitrary precision) integer object implementation */
  33.  
  34. /* XXX The functional organization of this file is terrible */
  35.  
  36. #include "allobjects.h"
  37. #include "longintrepr.h"
  38. #include "mymath.h"
  39. #include <assert.h>
  40. #include <ctype.h>
  41.  
  42. #define ABS(x) ((x) < 0 ? -(x) : (x))
  43.  
  44. /* Forward */
  45. #include "protos/longobject_protos.h"
  46. static longobject *long_normalize PROTO((longobject *));
  47. static longobject *mul1 PROTO((longobject *, wdigit));
  48. static longobject *muladd1 PROTO((longobject *, wdigit, wdigit));
  49. static longobject *divrem1 PROTO((longobject *, wdigit, digit *));
  50. static object *long_format PROTO((object *aa, int base));
  51.  
  52. static int ticker;    /* XXX Could be shared with ceval? */
  53.  
  54. #define SIGCHECK(block) \
  55.     if (--ticker < 0) { \
  56.         ticker = 100; \
  57.         if (sigcheck()) { block; } \
  58.     }
  59.  
  60. /* Normalize (remove leading zeros from) a long int object.
  61.    Doesn't attempt to free the storage--in most cases, due to the nature
  62.    of the algorithms used, this could save at most be one word anyway. */
  63.  
  64. static longobject *
  65. long_normalize(v)
  66.     register longobject *v;
  67. {
  68.     int j = ABS(v->ob_size);
  69.     register int i = j;
  70.     
  71.     while (i > 0 && v->ob_digit[i-1] == 0)
  72.         --i;
  73.     if (i != j)
  74.         v->ob_size = (v->ob_size < 0) ? -(i) : i;
  75.     return v;
  76. }
  77.  
  78. /* Allocate a new long int object with size digits.
  79.    Return NULL and set exception if we run out of memory. */
  80.  
  81. longobject *
  82. alloclongobject(size)
  83.     int size;
  84. {
  85.     return NEWVAROBJ(longobject, &Longtype, size);
  86. }
  87.  
  88. /* Create a new long int object from a C long int */
  89.  
  90. object *
  91. newlongobject(ival)
  92.     long ival;
  93. {
  94.     /* Assume a C long fits in at most 3 'digits' */
  95.     /* XXX On 64 bit machines this isn't true!!! */
  96.     longobject *v = alloclongobject(3);
  97.     if (v != NULL) {
  98.         if (ival < 0) {
  99.             ival = -ival;
  100.             v->ob_size = -(v->ob_size);
  101.         }
  102.         v->ob_digit[0] = ival & MASK;
  103.         v->ob_digit[1] = (ival >> SHIFT) & MASK;
  104.         v->ob_digit[2] = ((unsigned long)ival >> (2*SHIFT)) & MASK;
  105.         v = long_normalize(v);
  106.     }
  107.     return (object *)v;
  108. }
  109.  
  110. /* Create a new long int object from a C double */
  111.  
  112. object *
  113. #ifdef MPW
  114. dnewlongobject(double dval)
  115. #else
  116. dnewlongobject(dval)
  117.     double dval;
  118. #endif /* MPW */
  119. {
  120.     longobject *v;
  121.     double frac;
  122.     int i, ndig, expo, neg;
  123.     neg = 0;
  124.     if (dval < 0.0) {
  125.         neg = 1;
  126.         dval = -dval;
  127.     }
  128.     frac = frexp(dval, &expo); /* dval = frac*2**expo; 0.0 <= frac < 1.0 */
  129.  
  130.     if (expo <= 0)
  131.         return newlongobject(0L);
  132.     ndig = (expo-1) / SHIFT + 1; /* Number of 'digits' in result */
  133.     v = alloclongobject(ndig);
  134.     if (v == NULL)
  135.         return NULL;
  136.     frac = ldexp(frac, (expo-1) % SHIFT + 1);
  137.     for (i = ndig; --i >= 0; ) {
  138.         long bits = (long)frac;
  139.         v->ob_digit[i] = bits;
  140.         frac = frac - (double)bits;
  141.         frac = ldexp(frac, SHIFT);
  142.     }
  143.     if (neg)
  144.         v->ob_size = -(v->ob_size);
  145.     return (object *)v;
  146. }
  147.  
  148. /* Get a C long int from a long int object.
  149.    Returns -1 and sets an error condition if overflow occurs. */
  150.  
  151. long
  152. getlongvalue(vv)
  153.     object *vv;
  154. {
  155.     register longobject *v;
  156.     long x, prev;
  157.     int i, sign;
  158.     
  159.     if (vv == NULL || !is_longobject(vv)) {
  160.         err_badcall();
  161.         return -1;
  162.     }
  163.     v = (longobject *)vv;
  164.     i = v->ob_size;
  165.     sign = 1;
  166.     x = 0;
  167.     if (i < 0) {
  168.         sign = -1;
  169.         i = -(i);
  170.     }
  171.     while (--i >= 0) {
  172.         prev = x;
  173.         x = (x << SHIFT) + v->ob_digit[i];
  174.         if ((x >> SHIFT) != prev) {
  175.             err_setstr(OverflowError,
  176.                 "long int too long to convert");
  177.             return -1;
  178.         }
  179.     }
  180.     return x * sign;
  181. }
  182.  
  183. /* Get a C double from a long int object.  No overflow check. */
  184.  
  185. double
  186. dgetlongvalue(vv)
  187.     object *vv;
  188. {
  189.     register longobject *v;
  190.     double x;
  191.     double multiplier = (double) (1L << SHIFT);
  192.     int i, sign;
  193.     
  194.     if (vv == NULL || !is_longobject(vv)) {
  195.         err_badcall();
  196.         return -1;
  197.     }
  198.     v = (longobject *)vv;
  199.     i = v->ob_size;
  200.     sign = 1;
  201.     x = 0.0;
  202.     if (i < 0) {
  203.         sign = -1;
  204.         i = -(i);
  205.     }
  206.     while (--i >= 0) {
  207.         x = x*multiplier + (double)v->ob_digit[i];
  208.     }
  209.     return x * sign;
  210. }
  211.  
  212. /* Multiply by a single digit, ignoring the sign. */
  213.  
  214. static longobject *
  215. mul1(a, n)
  216.     longobject *a;
  217.     wdigit n;
  218. {
  219.     return muladd1(a, n, (digit)0);
  220. }
  221.  
  222. /* Multiply by a single digit and add a single digit, ignoring the sign. */
  223.  
  224. static longobject *
  225. muladd1(a, n, extra)
  226.     longobject *a;
  227.     wdigit n;
  228.     wdigit extra;
  229. {
  230.     int size_a = ABS(a->ob_size);
  231.     longobject *z = alloclongobject(size_a+1);
  232.     twodigits carry = extra;
  233.     int i;
  234.     
  235.     if (z == NULL)
  236.         return NULL;
  237.     for (i = 0; i < size_a; ++i) {
  238.         carry += (twodigits)a->ob_digit[i] * n;
  239.         z->ob_digit[i] = carry & MASK;
  240.         carry >>= SHIFT;
  241.     }
  242.     z->ob_digit[i] = carry;
  243.     return long_normalize(z);
  244. }
  245.  
  246. /* Divide a long integer by a digit, returning both the quotient
  247.    (as function result) and the remainder (through *prem).
  248.    The sign of a is ignored; n should not be zero. */
  249.  
  250. static longobject *
  251. divrem1(a, n, prem)
  252.     longobject *a;
  253.     wdigit n;
  254.     digit *prem;
  255. {
  256.     int size = ABS(a->ob_size);
  257.     longobject *z;
  258.     int i;
  259.     twodigits rem = 0;
  260.     
  261.     assert(n > 0 && n <= MASK);
  262.     z = alloclongobject(size);
  263.     if (z == NULL)
  264.         return NULL;
  265.     for (i = size; --i >= 0; ) {
  266.         rem = (rem << SHIFT) + a->ob_digit[i];
  267.         z->ob_digit[i] = rem/n;
  268.         rem %= n;
  269.     }
  270.     *prem = rem;
  271.     return long_normalize(z);
  272. }
  273.  
  274. /* Convert a long int object to a string, using a given conversion base.
  275.    Return a string object.
  276.    If base is 8 or 16, add the proper prefix '0' or '0x'.
  277.    External linkage: used in bltinmodule.c by hex() and oct(). */
  278.  
  279. static object *
  280. long_format(aa, base)
  281.     object *aa;
  282.     int base;
  283. {
  284.     register longobject *a = (longobject *)aa;
  285.     stringobject *str;
  286.     int i;
  287.     int size_a = ABS(a->ob_size);
  288.     char *p;
  289.     int bits;
  290.     char sign = '\0';
  291.  
  292.     if (a == NULL || !is_longobject(a)) {
  293.         err_badcall();
  294.         return NULL;
  295.     }
  296.     assert(base >= 2 && base <= 36);
  297.     
  298.     /* Compute a rough upper bound for the length of the string */
  299.     i = base;
  300.     bits = 0;
  301.     while (i > 1) {
  302.         ++bits;
  303.         i >>= 1;
  304.     }
  305.     i = 6 + (size_a*SHIFT + bits-1) / bits;
  306.     str = (stringobject *) newsizedstringobject((char *)0, i);
  307.     if (str == NULL)
  308.         return NULL;
  309.     p = GETSTRINGVALUE(str) + i;
  310.     *p = '\0';
  311.     *--p = 'L';
  312.     if (a->ob_size < 0)
  313.         sign = '-';
  314.     
  315.     INCREF(a);
  316.     do {
  317.         digit rem;
  318.         longobject *temp = divrem1(a, (digit)base, &rem);
  319.         if (temp == NULL) {
  320.             DECREF(a);
  321.             DECREF(str);
  322.             return NULL;
  323.         }
  324.         if (rem < 10)
  325.             rem += '0';
  326.         else
  327.             rem += 'A'-10;
  328.         assert(p > GETSTRINGVALUE(str));
  329.         *--p = rem;
  330.         DECREF(a);
  331.         a = temp;
  332.         SIGCHECK({
  333.             DECREF(a);
  334.             DECREF(str);
  335.             return NULL;
  336.         })
  337.     } while (ABS(a->ob_size) != 0);
  338.     DECREF(a);
  339.     if (base == 8) {
  340.         if (size_a != 0)
  341.             *--p = '0';
  342.     }
  343.     else if (base == 16) {
  344.         *--p = 'x';
  345.         *--p = '0';
  346.     }
  347.     else if (base != 10) {
  348.         *--p = '#';
  349.         *--p = '0' + base%10;
  350.         if (base > 10)
  351.             *--p = '0' + base/10;
  352.     }
  353.     if (sign)
  354.         *--p = sign;
  355.     if (p != GETSTRINGVALUE(str)) {
  356.         char *q = GETSTRINGVALUE(str);
  357.         assert(p > q);
  358.         do {
  359.         } while ((*q++ = *p++) != '\0');
  360.         q--;
  361.         resizestring((object **)&str, (int) (q - GETSTRINGVALUE(str)));
  362.     }
  363.     return (object *)str;
  364. }
  365.  
  366. #if 0
  367. /* Convert a string to a long int object, in a given base.
  368.    Base zero implies a default depending on the number.
  369.    External linkage: used in compile.c and stropmodule.c. */
  370.  
  371. object *
  372. long_scan(str, base)
  373.     char *str;
  374.     int base;
  375. {
  376.     return long_escan(str, (char **)NULL, base);
  377. }
  378. #endif
  379.  
  380. object *
  381. long_escan(str, pend, base)
  382.     char *str;
  383.     char **pend;
  384.     int base;
  385. {
  386.     int sign = 1;
  387.     longobject *z;
  388.     
  389.     if (base != 0 && base < 2 || base > 36) {
  390.         err_setstr(ValueError, "invalid base for long literal");
  391.         return NULL;
  392.     }
  393.     while (*str != '\0' && isspace(Py_CHARMASK(*str)))
  394.         str++;
  395.     if (*str == '+')
  396.         ++str;
  397.     else if (*str == '-') {
  398.         ++str;
  399.         sign = -1;
  400.     }
  401.     while (*str != '\0' && isspace(Py_CHARMASK(*str)))
  402.         str++;
  403.     if (base == 0) {
  404.         if (str[0] != '0')
  405.             base = 10;
  406.         else if (str[1] == 'x' || str[1] == 'X')
  407.             base = 16;
  408.         else
  409.             base = 8;
  410.     }
  411.     if (base == 16 && str[0] == '0' && (str[1] == 'x' || str[1] == 'X'))
  412.         str += 2;
  413.     z = alloclongobject(0);
  414.     for ( ; z != NULL; ++str) {
  415.         int k = -1;
  416.         longobject *temp;
  417.         
  418.         if (*str <= '9')
  419.             k = *str - '0';
  420.         else if (*str >= 'a')
  421.             k = *str - 'a' + 10;
  422.         else if (*str >= 'A')
  423.             k = *str - 'A' + 10;
  424.         if (k < 0 || k >= base)
  425.             break;
  426.         temp = muladd1(z, (digit)base, (digit)k);
  427.         DECREF(z);
  428.         z = temp;
  429.     }
  430.     if (sign < 0 && z != NULL && z->ob_size != 0)
  431.         z->ob_size = -(z->ob_size);
  432.     if (pend)
  433.         *pend = str;
  434.     return (object *) z;
  435. }
  436.  
  437. static longobject *x_divrem PROTO((longobject *, longobject *, longobject **));
  438. static object *long_pos PROTO((longobject *));
  439. static long_divrem PROTO((longobject *, longobject *,
  440.     longobject **, longobject **));
  441.  
  442. /* Long division with remainder, top-level routine */
  443.  
  444. static int
  445. long_divrem(a, b, pdiv, prem)
  446.     longobject *a, *b;
  447.     longobject **pdiv;
  448.     longobject **prem;
  449. {
  450.     int size_a = ABS(a->ob_size), size_b = ABS(b->ob_size);
  451.     longobject *z;
  452.     
  453.     if (size_b == 0) {
  454.         err_setstr(ZeroDivisionError, "long division or modulo");
  455.         return -1;
  456.     }
  457.     if (size_a < size_b ||
  458.             size_a == size_b &&
  459.             a->ob_digit[size_a-1] < b->ob_digit[size_b-1]) {
  460.         /* |a| < |b|. */
  461.         *pdiv = alloclongobject(0);
  462.         INCREF(a);
  463.         *prem = (longobject *) a;
  464.         return 0;
  465.     }
  466.     if (size_b == 1) {
  467.         digit rem = 0;
  468.         z = divrem1(a, b->ob_digit[0], &rem);
  469.         if (z == NULL)
  470.             return -1;
  471.         *prem = (longobject *) newlongobject((long)rem);
  472.     }
  473.     else {
  474.         z = x_divrem(a, b, prem);
  475.         if (z == NULL)
  476.             return -1;
  477.     }
  478.     /* Set the signs.
  479.        The quotient z has the sign of a*b;
  480.        the remainder r has the sign of a,
  481.        so a = b*z + r. */
  482.     if ((a->ob_size < 0) != (b->ob_size < 0))
  483.         z->ob_size = -(z->ob_size);
  484.     if (a->ob_size < 0 && (*prem)->ob_size != 0)
  485.         (*prem)->ob_size = -((*prem)->ob_size);
  486.     *pdiv = z;
  487.     return 0;
  488. }
  489.  
  490. /* Unsigned long division with remainder -- the algorithm */
  491.  
  492. static longobject *
  493. x_divrem(v1, w1, prem)
  494.     longobject *v1, *w1;
  495.     longobject **prem;
  496. {
  497.     int size_v = ABS(v1->ob_size), size_w = ABS(w1->ob_size);
  498.     digit d = (twodigits)BASE / (w1->ob_digit[size_w-1] + 1);
  499.     longobject *v = mul1(v1, d);
  500.     longobject *w = mul1(w1, d);
  501.     longobject *a;
  502.     int j, k;
  503.     
  504.     if (v == NULL || w == NULL) {
  505.         XDECREF(v);
  506.         XDECREF(w);
  507.         return NULL;
  508.     }
  509.     
  510.     assert(size_v >= size_w && size_w > 1); /* Assert checks by div() */
  511.     assert(v->ob_refcnt == 1); /* Since v will be used as accumulator! */
  512.     assert(size_w == ABS(w->ob_size)); /* That's how d was calculated */
  513.     
  514.     size_v = ABS(v->ob_size);
  515.     a = alloclongobject(size_v - size_w + 1);
  516.     
  517.     for (j = size_v, k = a->ob_size-1; a != NULL && k >= 0; --j, --k) {
  518.         digit vj = (j >= size_v) ? 0 : v->ob_digit[j];
  519.         twodigits q;
  520.         stwodigits carry = 0;
  521.         int i;
  522.         
  523.         SIGCHECK({
  524.             DECREF(a);
  525.             a = NULL;
  526.             break;
  527.         })
  528.         if (vj == w->ob_digit[size_w-1])
  529.             q = MASK;
  530.         else
  531.             q = (((twodigits)vj << SHIFT) + v->ob_digit[j-1]) /
  532.                 w->ob_digit[size_w-1];
  533.         
  534.         while (w->ob_digit[size_w-2]*q >
  535.                 ((
  536.                     ((twodigits)vj << SHIFT)
  537.                     + v->ob_digit[j-1]
  538.                     - q*w->ob_digit[size_w-1]
  539.                                 ) << SHIFT)
  540.                 + v->ob_digit[j-2])
  541.             --q;
  542.         
  543.         for (i = 0; i < size_w && i+k < size_v; ++i) {
  544.             twodigits z = w->ob_digit[i] * q;
  545.             digit zz = z >> SHIFT;
  546.             carry += v->ob_digit[i+k] - z + ((twodigits)zz << SHIFT);
  547.             v->ob_digit[i+k] = carry & MASK;
  548.             carry = (carry >> SHIFT) - zz;
  549.         }
  550.         
  551.         if (i+k < size_v) {
  552.             carry += v->ob_digit[i+k];
  553.             v->ob_digit[i+k] = 0;
  554.         }
  555.         
  556.         if (carry == 0)
  557.             a->ob_digit[k] = q;
  558.         else {
  559.             assert(carry == -1);
  560.             a->ob_digit[k] = q-1;
  561.             carry = 0;
  562.             for (i = 0; i < size_w && i+k < size_v; ++i) {
  563.                 carry += v->ob_digit[i+k] + w->ob_digit[i];
  564.                 v->ob_digit[i+k] = carry & MASK;
  565.                 carry >>= SHIFT;
  566.             }
  567.         }
  568.     } /* for j, k */
  569.     
  570.     if (a == NULL)
  571.         *prem = NULL;
  572.     else {
  573.         a = long_normalize(a);
  574.         *prem = divrem1(v, d, &d);
  575.         /* d receives the (unused) remainder */
  576.         if (*prem == NULL) {
  577.             DECREF(a);
  578.             a = NULL;
  579.         }
  580.     }
  581.     DECREF(v);
  582.     DECREF(w);
  583.     return a;
  584. }
  585.  
  586. /* Methods */
  587.  
  588. /* Forward */
  589. static void long_dealloc PROTO((object *));
  590. static object *long_repr PROTO((object *));
  591. static int long_compare PROTO((longobject *, longobject *));
  592. static long long_hash PROTO((longobject *));
  593.  
  594. static object *long_add PROTO((longobject *, longobject *));
  595. static object *long_sub PROTO((longobject *, longobject *));
  596. static object *long_mul PROTO((longobject *, longobject *));
  597. static object *long_div PROTO((longobject *, longobject *));
  598. static object *long_mod PROTO((longobject *, longobject *));
  599. static object *long_divmod PROTO((longobject *, longobject *));
  600. static object *long_pow PROTO((longobject *, longobject *, longobject *));
  601. static object *long_neg PROTO((longobject *));
  602. static object *long_pos PROTO((longobject *));
  603. static object *long_abs PROTO((longobject *));
  604. static int long_nonzero PROTO((longobject *));
  605. static object *long_invert PROTO((longobject *));
  606. static object *long_lshift PROTO((longobject *, longobject *));
  607. static object *long_rshift PROTO((longobject *, longobject *));
  608. static object *long_and PROTO((longobject *, longobject *));
  609. static object *long_xor PROTO((longobject *, longobject *));
  610. static object *long_or PROTO((longobject *, longobject *));
  611.  
  612. static void
  613. long_dealloc(v)
  614.     object *v;
  615. {
  616.     DEL(v);
  617. }
  618.  
  619. static object *
  620. long_repr(v)
  621.     object *v;
  622. {
  623.     return long_format(v, 10);
  624. }
  625.  
  626. static int
  627. long_compare(a, b)
  628.     longobject *a, *b;
  629. {
  630.     int sign;
  631.     
  632.     if (a->ob_size != b->ob_size) {
  633.         if (ABS(a->ob_size) == 0 && ABS(b->ob_size) == 0)
  634.             sign = 0;
  635.         else
  636.             sign = a->ob_size - b->ob_size;
  637.     }
  638.     else {
  639.         int i = ABS(a->ob_size);
  640.         while (--i >= 0 && a->ob_digit[i] == b->ob_digit[i])
  641.             ;
  642.         if (i < 0)
  643.             sign = 0;
  644.         else {
  645.             sign = (int)a->ob_digit[i] - (int)b->ob_digit[i];
  646.             if (a->ob_size < 0)
  647.                 sign = -sign;
  648.         }
  649.     }
  650.     return sign < 0 ? -1 : sign > 0 ? 1 : 0;
  651. }
  652.  
  653. static long
  654. long_hash(v)
  655.     longobject *v;
  656. {
  657.     long x;
  658.     int i, sign;
  659.  
  660.     /* This is designed so that Python ints and longs with the
  661.        same value hash to the same value, otherwise comparisons
  662.        of mapping keys will turn out weird */
  663.     i = v->ob_size;
  664.     sign = 1;
  665.     x = 0;
  666.     if (i < 0) {
  667.         sign = -1;
  668.         i = -(i);
  669.     }
  670.     while (--i >= 0) {
  671.         /* Force a 32-bit circular shift */
  672.         x = ((x << SHIFT) & ~MASK) | ((x >> (32-SHIFT)) & MASK);
  673.         x += v->ob_digit[i];
  674.     }
  675.     x = x * sign;
  676.     if (x == -1)
  677.         x = -2;
  678.     return x;
  679. }
  680.  
  681.  
  682. /* Add the absolute values of two long integers. */
  683.  
  684. static longobject *x_add PROTO((longobject *, longobject *));
  685. static longobject *
  686. x_add(a, b)
  687.     longobject *a, *b;
  688. {
  689.     int size_a = ABS(a->ob_size), size_b = ABS(b->ob_size);
  690.     longobject *z;
  691.     int i;
  692.     digit carry = 0;
  693.     
  694.     /* Ensure a is the larger of the two: */
  695.     if (size_a < size_b) {
  696.         { longobject *temp = a; a = b; b = temp; }
  697.         { int size_temp = size_a; size_a = size_b; size_b = size_temp; }
  698.     }
  699.     z = alloclongobject(size_a+1);
  700.     if (z == NULL)
  701.         return NULL;
  702.     for (i = 0; i < size_b; ++i) {
  703.         carry += a->ob_digit[i] + b->ob_digit[i];
  704.         z->ob_digit[i] = carry & MASK;
  705.         /* The following assumes unsigned shifts don't
  706.            propagate the sign bit. */
  707.         carry >>= SHIFT;
  708.     }
  709.     for (; i < size_a; ++i) {
  710.         carry += a->ob_digit[i];
  711.         z->ob_digit[i] = carry & MASK;
  712.         carry >>= SHIFT;
  713.     }
  714.     z->ob_digit[i] = carry;
  715.     return long_normalize(z);
  716. }
  717.  
  718. /* Subtract the absolute values of two integers. */
  719.  
  720. static longobject *x_sub PROTO((longobject *, longobject *));
  721. static longobject *
  722. x_sub(a, b)
  723.     longobject *a, *b;
  724. {
  725.     int size_a = ABS(a->ob_size), size_b = ABS(b->ob_size);
  726.     longobject *z;
  727.     int i;
  728.     int sign = 1;
  729.     digit borrow = 0;
  730.     
  731.     /* Ensure a is the larger of the two: */
  732.     if (size_a < size_b) {
  733.         sign = -1;
  734.         { longobject *temp = a; a = b; b = temp; }
  735.         { int size_temp = size_a; size_a = size_b; size_b = size_temp; }
  736.     }
  737.     else if (size_a == size_b) {
  738.         /* Find highest digit where a and b differ: */
  739.         i = size_a;
  740.         while (--i >= 0 && a->ob_digit[i] == b->ob_digit[i])
  741.             ;
  742.         if (i < 0)
  743.             return alloclongobject(0);
  744.         if (a->ob_digit[i] < b->ob_digit[i]) {
  745.             sign = -1;
  746.             { longobject *temp = a; a = b; b = temp; }
  747.         }
  748.         size_a = size_b = i+1;
  749.     }
  750.     z = alloclongobject(size_a);
  751.     if (z == NULL)
  752.         return NULL;
  753.     for (i = 0; i < size_b; ++i) {
  754.         /* The following assumes unsigned arithmetic
  755.            works module 2**N for some N>SHIFT. */
  756.         borrow = a->ob_digit[i] - b->ob_digit[i] - borrow;
  757.         z->ob_digit[i] = borrow & MASK;
  758.         borrow >>= SHIFT;
  759.         borrow &= 1; /* Keep only one sign bit */
  760.     }
  761.     for (; i < size_a; ++i) {
  762.         borrow = a->ob_digit[i] - borrow;
  763.         z->ob_digit[i] = borrow & MASK;
  764.         borrow >>= SHIFT;
  765.     }
  766.     assert(borrow == 0);
  767.     if (sign < 0)
  768.         z->ob_size = -(z->ob_size);
  769.     return long_normalize(z);
  770. }
  771.  
  772. static object *
  773. long_add(a, b)
  774.     longobject *a;
  775.     longobject *b;
  776. {
  777.     longobject *z;
  778.     
  779.     if (a->ob_size < 0) {
  780.         if (b->ob_size < 0) {
  781.             z = x_add(a, b);
  782.             if (z != NULL && z->ob_size != 0)
  783.                 z->ob_size = -(z->ob_size);
  784.         }
  785.         else
  786.             z = x_sub(b, a);
  787.     }
  788.     else {
  789.         if (b->ob_size < 0)
  790.             z = x_sub(a, b);
  791.         else
  792.             z = x_add(a, b);
  793.     }
  794.     return (object *)z;
  795. }
  796.  
  797. static object *
  798. long_sub(a, b)
  799.     longobject *a;
  800.     longobject *b;
  801. {
  802.     longobject *z;
  803.     
  804.     if (a->ob_size < 0) {
  805.         if (b->ob_size < 0)
  806.             z = x_sub(a, b);
  807.         else
  808.             z = x_add(a, b);
  809.         if (z != NULL && z->ob_size != 0)
  810.             z->ob_size = -(z->ob_size);
  811.     }
  812.     else {
  813.         if (b->ob_size < 0)
  814.             z = x_add(a, b);
  815.         else
  816.             z = x_sub(a, b);
  817.     }
  818.     return (object *)z;
  819. }
  820.  
  821. static object *
  822. long_mul(a, b)
  823.     longobject *a;
  824.     longobject *b;
  825. {
  826.     int size_a;
  827.     int size_b;
  828.     longobject *z;
  829.     int i;
  830.     
  831.     size_a = ABS(a->ob_size);
  832.     size_b = ABS(b->ob_size);
  833.     z = alloclongobject(size_a + size_b);
  834.     if (z == NULL)
  835.         return NULL;
  836.     for (i = 0; i < z->ob_size; ++i)
  837.         z->ob_digit[i] = 0;
  838.     for (i = 0; i < size_a; ++i) {
  839.         twodigits carry = 0;
  840.         twodigits f = a->ob_digit[i];
  841.         int j;
  842.         
  843.         SIGCHECK({
  844.             DECREF(z);
  845.             return NULL;
  846.         })
  847.         for (j = 0; j < size_b; ++j) {
  848.             carry += z->ob_digit[i+j] + b->ob_digit[j] * f;
  849.             z->ob_digit[i+j] = carry & MASK;
  850.             carry >>= SHIFT;
  851.         }
  852.         for (; carry != 0; ++j) {
  853.             assert(i+j < z->ob_size);
  854.             carry += z->ob_digit[i+j];
  855.             z->ob_digit[i+j] = carry & MASK;
  856.             carry >>= SHIFT;
  857.         }
  858.     }
  859.     if (a->ob_size < 0)
  860.         z->ob_size = -(z->ob_size);
  861.     if (b->ob_size < 0)
  862.         z->ob_size = -(z->ob_size);
  863.     return (object *) long_normalize(z);
  864. }
  865.  
  866. /* The / and % operators are now defined in terms of divmod().
  867.    The expression a mod b has the value a - b*floor(a/b).
  868.    The long_divrem function gives the remainder after division of
  869.    |a| by |b|, with the sign of a.  This is also expressed
  870.    as a - b*trunc(a/b), if trunc truncates towards zero.
  871.    Some examples:
  872.         a     b    a rem b        a mod b
  873.         13     10     3         3
  874.        -13     10    -3         7
  875.         13    -10     3        -7
  876.        -13    -10    -3        -3
  877.    So, to get from rem to mod, we have to add b if a and b
  878.    have different signs.  We then subtract one from the 'div'
  879.    part of the outcome to keep the invariant intact. */
  880.  
  881. static int l_divmod PROTO((longobject *, longobject *,
  882.     longobject **, longobject **));
  883. static int
  884. l_divmod(v, w, pdiv, pmod)
  885.     longobject *v;
  886.     longobject *w;
  887.     longobject **pdiv;
  888.     longobject **pmod;
  889. {
  890.     longobject *div, *mod;
  891.     
  892.     if (long_divrem(v, w, &div, &mod) < 0)
  893.         return -1;
  894.     if (mod->ob_size < 0 && w->ob_size > 0 ||
  895.                 mod->ob_size > 0 && w->ob_size < 0) {
  896.         longobject *temp;
  897.         longobject *one;
  898.         temp = (longobject *) long_add(mod, w);
  899.         DECREF(mod);
  900.         mod = temp;
  901.         if (mod == NULL) {
  902.             DECREF(div);
  903.             return -1;
  904.         }
  905.         one = (longobject *) newlongobject(1L);
  906.         if (one == NULL ||
  907.             (temp = (longobject *) long_sub(div, one)) == NULL) {
  908.             DECREF(mod);
  909.             DECREF(div);
  910.             XDECREF(one);
  911.             return -1;
  912.         }
  913.         DECREF(one);
  914.         DECREF(div);
  915.         div = temp;
  916.     }
  917.     *pdiv = div;
  918.     *pmod = mod;
  919.     return 0;
  920. }
  921.  
  922. static object *
  923. long_div(v, w)
  924.     longobject *v;
  925.     longobject *w;
  926. {
  927.     longobject *div, *mod;
  928.     if (l_divmod(v, w, &div, &mod) < 0)
  929.         return NULL;
  930.     DECREF(mod);
  931.     return (object *)div;
  932. }
  933.  
  934. static object *
  935. long_mod(v, w)
  936.     longobject *v;
  937.     longobject *w;
  938. {
  939.     longobject *div, *mod;
  940.     if (l_divmod(v, w, &div, &mod) < 0)
  941.         return NULL;
  942.     DECREF(div);
  943.     return (object *)mod;
  944. }
  945.  
  946. static object *
  947. long_divmod(v, w)
  948.     longobject *v;
  949.     longobject *w;
  950. {
  951.     object *z;
  952.     longobject *div, *mod;
  953.     if (l_divmod(v, w, &div, &mod) < 0)
  954.         return NULL;
  955.     z = newtupleobject(2);
  956.     if (z != NULL) {
  957.         settupleitem(z, 0, (object *) div);
  958.         settupleitem(z, 1, (object *) mod);
  959.     }
  960.     else {
  961.         DECREF(div);
  962.         DECREF(mod);
  963.     }
  964.     return z;
  965. }
  966.  
  967. static object *
  968. long_pow(a, b, c)
  969.     longobject *a;
  970.     longobject *b;
  971.     longobject *c;
  972. {
  973.     longobject *z, *div, *mod;
  974.     int size_b, i;
  975.     
  976.     size_b = b->ob_size;
  977.     if (size_b < 0) {
  978.         err_setstr(ValueError, "long integer to the negative power");
  979.         return NULL;
  980.     }
  981.     z = (longobject *)newlongobject(1L);
  982.     INCREF(a);
  983.     for (i = 0; i < size_b; ++i) {
  984.         digit bi = b->ob_digit[i];
  985.         int j;
  986.     
  987.         for (j = 0; j < SHIFT; ++j) {
  988.             longobject *temp;
  989.         
  990.             if (bi & 1) {
  991.                 temp = (longobject *)long_mul(z, a);
  992.                 DECREF(z);
  993.                  if ((object*)c!=None && temp!=NULL) {
  994.                      l_divmod(temp, c, &div, &mod);
  995.                      XDECREF(div);
  996.                      DECREF(temp);
  997.                      temp = mod;
  998.                 }
  999.                  z = temp;
  1000.                 if (z == NULL)
  1001.                     break;
  1002.             }
  1003.             bi >>= 1;
  1004.             if (bi == 0 && i+1 == size_b)
  1005.                 break;
  1006.             temp = (longobject *)long_mul(a, a);
  1007.             DECREF(a);
  1008.              if ((object*)c!=None && temp!=NULL) {
  1009.                  l_divmod(temp, c, &div, &mod);
  1010.                  XDECREF(div);
  1011.                  DECREF(temp);
  1012.                  temp = mod;
  1013.             }
  1014.             a = temp;
  1015.             if (a == NULL) {
  1016.                 DECREF(z);
  1017.                 z = NULL;
  1018.                 break;
  1019.             }
  1020.         }
  1021.         if (a == NULL || z == NULL)
  1022.             break;
  1023.     }
  1024.     XDECREF(a);
  1025.     if ((object*)c!=None && z!=NULL) {
  1026.             l_divmod(z, c, &div, &mod);
  1027.             XDECREF(div);
  1028.             DECREF(z);
  1029.             z=mod;
  1030.     }
  1031.     return (object *)z;
  1032. }
  1033.  
  1034. static object *
  1035. long_invert(v)
  1036.     longobject *v;
  1037. {
  1038.     /* Implement ~x as -(x+1) */
  1039.     longobject *x;
  1040.     longobject *w;
  1041.     w = (longobject *)newlongobject(1L);
  1042.     if (w == NULL)
  1043.         return NULL;
  1044.     x = (longobject *) long_add(v, w);
  1045.     DECREF(w);
  1046.     if (x == NULL)
  1047.         return NULL;
  1048.     if (x->ob_size != 0)
  1049.         x->ob_size = -(x->ob_size);
  1050.     return (object *)x;
  1051. }
  1052.  
  1053. static object *
  1054. long_pos(v)
  1055.     longobject *v;
  1056. {
  1057.     INCREF(v);
  1058.     return (object *)v;
  1059. }
  1060.  
  1061. static object *
  1062. long_neg(v)
  1063.     longobject *v;
  1064. {
  1065.     longobject *z;
  1066.     int i, n;
  1067.     n = ABS(v->ob_size);
  1068.     if (n == 0) {
  1069.         /* -0 == 0 */
  1070.         INCREF(v);
  1071.         return (object *) v;
  1072.     }
  1073.     z = alloclongobject(ABS(n));
  1074.     if (z == NULL)
  1075.         return NULL;
  1076.     for (i = 0; i < n; i++)
  1077.         z->ob_digit[i] = v->ob_digit[i];
  1078.     z->ob_size = -(v->ob_size);
  1079.     return (object *)z;
  1080. }
  1081.  
  1082. static object *
  1083. long_abs(v)
  1084.     longobject *v;
  1085. {
  1086.     if (v->ob_size < 0)
  1087.         return long_neg(v);
  1088.     else {
  1089.         INCREF(v);
  1090.         return (object *)v;
  1091.     }
  1092. }
  1093.  
  1094. static int
  1095. long_nonzero(v)
  1096.     longobject *v;
  1097. {
  1098.     return ABS(v->ob_size) != 0;
  1099. }
  1100.  
  1101. static object *
  1102. long_rshift(a, b)
  1103.     longobject *a;
  1104.     longobject *b;
  1105. {
  1106.     longobject *z;
  1107.     long shiftby;
  1108.     int newsize, wordshift, loshift, hishift, i, j;
  1109.     digit lomask, himask;
  1110.     
  1111.     if (a->ob_size < 0) {
  1112.         /* Right shifting negative numbers is harder */
  1113.         longobject *a1, *a2, *a3;
  1114.         a1 = (longobject *) long_invert(a);
  1115.         if (a1 == NULL) return NULL;
  1116.         a2 = (longobject *) long_rshift(a1, b);
  1117.         DECREF(a1);
  1118.         if (a2 == NULL) return NULL;
  1119.         a3 = (longobject *) long_invert(a2);
  1120.         DECREF(a2);
  1121.         return (object *) a3;
  1122.     }
  1123.     
  1124.     shiftby = getlongvalue((object *)b);
  1125.     if (shiftby == -1L && err_occurred())
  1126.         return NULL;
  1127.     if (shiftby < 0) {
  1128.         err_setstr(ValueError, "negative shift count");
  1129.         return NULL;
  1130.     }
  1131.     wordshift = shiftby / SHIFT;
  1132.     newsize = ABS(a->ob_size) - wordshift;
  1133.     if (newsize <= 0) {
  1134.         z = alloclongobject(0);
  1135.         return (object *)z;
  1136.     }
  1137.     loshift = shiftby % SHIFT;
  1138.     hishift = SHIFT - loshift;
  1139.     lomask = ((digit)1 << hishift) - 1;
  1140.     himask = MASK ^ lomask;
  1141.     z = alloclongobject(newsize);
  1142.     if (z == NULL)
  1143.         return NULL;
  1144.     if (a->ob_size < 0)
  1145.         z->ob_size = -(z->ob_size);
  1146.     for (i = 0, j = wordshift; i < newsize; i++, j++) {
  1147.         z->ob_digit[i] = (a->ob_digit[j] >> loshift) & lomask;
  1148.         if (i+1 < newsize)
  1149.             z->ob_digit[i] |=
  1150.               (a->ob_digit[j+1] << hishift) & himask;
  1151.     }
  1152.     return (object *) long_normalize(z);
  1153. }
  1154.  
  1155. static object *
  1156. long_lshift(a, b)
  1157.     longobject *a;
  1158.     longobject *b;
  1159. {
  1160.     longobject *z;
  1161.     long shiftby;
  1162.     int newsize, wordshift, loshift, hishift, i, j;
  1163.     digit lomask, himask;
  1164.     
  1165.     shiftby = getlongvalue((object *)b);
  1166.     if (shiftby == -1L && err_occurred())
  1167.         return NULL;
  1168.     if (shiftby < 0) {
  1169.         err_setstr(ValueError, "negative shift count");
  1170.         return NULL;
  1171.     }
  1172.     if (shiftby > MASK) {
  1173.         err_setstr(ValueError, "outrageous left shift count");
  1174.         return NULL;
  1175.     }
  1176.     if (shiftby % SHIFT == 0) {
  1177.         wordshift = shiftby / SHIFT;
  1178.         loshift = 0;
  1179.         hishift = SHIFT;
  1180.         newsize = ABS(a->ob_size) + wordshift;
  1181.         lomask = MASK;
  1182.         himask = 0;
  1183.     }
  1184.     else {
  1185.         wordshift = shiftby / SHIFT + 1;
  1186.         loshift = SHIFT - shiftby%SHIFT;
  1187.         hishift = shiftby % SHIFT;
  1188.         newsize = ABS(a->ob_size) + wordshift;
  1189.         lomask = ((digit)1 << hishift) - 1;
  1190.         himask = MASK ^ lomask;
  1191.     }
  1192.     z = alloclongobject(newsize);
  1193.     if (z == NULL)
  1194.         return NULL;
  1195.     if (a->ob_size < 0)
  1196.         z->ob_size = -(z->ob_size);
  1197.     for (i = 0; i < wordshift; i++)
  1198.         z->ob_digit[i] = 0;
  1199.     for (i = wordshift, j = 0; i < newsize; i++, j++) {
  1200.         if (i > 0)
  1201.             z->ob_digit[i-1] |=
  1202.                 (a->ob_digit[j] << hishift) & himask;
  1203.         z->ob_digit[i] =
  1204.             (a->ob_digit[j] >> loshift) & lomask;
  1205.     }
  1206.     return (object *) long_normalize(z);
  1207. }
  1208.  
  1209.  
  1210. /* Bitwise and/xor/or operations */
  1211.  
  1212. #define MAX(x, y) ((x) < (y) ? (y) : (x))
  1213. #define MIN(x, y) ((x) > (y) ? (y) : (x))
  1214.  
  1215. static object *long_bitwise PROTO((longobject *, int, longobject *));
  1216. static object *
  1217. long_bitwise(a, op, b)
  1218.     longobject *a;
  1219.     int op; /* '&', '|', '^' */
  1220.     longobject *b;
  1221. {
  1222.     digit maska, maskb; /* 0 or MASK */
  1223.     int negz;
  1224.     int size_a, size_b, size_z;
  1225.     longobject *z;
  1226.     int i;
  1227.     digit diga, digb;
  1228.     object *v;
  1229.     
  1230.     if (a->ob_size < 0) {
  1231.         a = (longobject *) long_invert(a);
  1232.         maska = MASK;
  1233.     }
  1234.     else {
  1235.         INCREF(a);
  1236.         maska = 0;
  1237.     }
  1238.     if (b->ob_size < 0) {
  1239.         b = (longobject *) long_invert(b);
  1240.         maskb = MASK;
  1241.     }
  1242.     else {
  1243.         INCREF(b);
  1244.         maskb = 0;
  1245.     }
  1246.     
  1247.     size_a = a->ob_size;
  1248.     size_b = b->ob_size;
  1249.     size_z = MAX(size_a, size_b);
  1250.     z = alloclongobject(size_z);
  1251.     if (a == NULL || b == NULL || z == NULL) {
  1252.         XDECREF(a);
  1253.         XDECREF(b);
  1254.         XDECREF(z);
  1255.         return NULL;
  1256.     }
  1257.     
  1258.     negz = 0;
  1259.     switch (op) {
  1260.     case '^':
  1261.         if (maska != maskb) {
  1262.             maska ^= MASK;
  1263.             negz = -1;
  1264.         }
  1265.         break;
  1266.     case '&':
  1267.         if (maska && maskb) {
  1268.             op = '|';
  1269.             maska ^= MASK;
  1270.             maskb ^= MASK;
  1271.             negz = -1;
  1272.         }
  1273.         break;
  1274.     case '|':
  1275.         if (maska || maskb) {
  1276.             op = '&';
  1277.             maska ^= MASK;
  1278.             maskb ^= MASK;
  1279.             negz = -1;
  1280.         }
  1281.         break;
  1282.     }
  1283.     
  1284.     for (i = 0; i < size_z; ++i) {
  1285.         diga = (i < size_a ? a->ob_digit[i] : 0) ^ maska;
  1286.         digb = (i < size_b ? b->ob_digit[i] : 0) ^ maskb;
  1287.         switch (op) {
  1288.         case '&': z->ob_digit[i] = diga & digb; break;
  1289.         case '|': z->ob_digit[i] = diga | digb; break;
  1290.         case '^': z->ob_digit[i] = diga ^ digb; break;
  1291.         }
  1292.     }
  1293.     
  1294.     DECREF(a);
  1295.     DECREF(b);
  1296.     z = long_normalize(z);
  1297.     if (negz == 0)
  1298.         return (object *) z;
  1299.     v = long_invert(z);
  1300.     DECREF(z);
  1301.     return v;
  1302. }
  1303.  
  1304. static object *
  1305. long_and(a, b)
  1306.     longobject *a;
  1307.     longobject *b;
  1308. {
  1309.     return long_bitwise(a, '&', b);
  1310. }
  1311.  
  1312. static object *
  1313. long_xor(a, b)
  1314.     longobject *a;
  1315.     longobject *b;
  1316. {
  1317.     return long_bitwise(a, '^', b);
  1318. }
  1319.  
  1320. static object *
  1321. long_or(a, b)
  1322.     longobject *a;
  1323.     longobject *b;
  1324. {
  1325.     return long_bitwise(a, '|', b);
  1326. }
  1327.  
  1328. static int
  1329. long_coerce(pv, pw)
  1330.     object **pv;
  1331.     object **pw;
  1332. {
  1333.     if (is_intobject(*pw)) {
  1334.         *pw = newlongobject(getintvalue(*pw));
  1335.         INCREF(*pv);
  1336.         return 0;
  1337.     }
  1338.     return 1; /* Can't do it */
  1339. }
  1340.  
  1341. static object *
  1342. long_int(v)
  1343.     object *v;
  1344. {
  1345.     long x;
  1346.     x = getlongvalue(v);
  1347.     if (err_occurred())
  1348.         return NULL;
  1349.     return newintobject(x);
  1350. }
  1351.  
  1352. static object *
  1353. long_long(v)
  1354.     object *v;
  1355. {
  1356.     INCREF(v);
  1357.     return v;
  1358. }
  1359.  
  1360. static object *
  1361. long_float(v)
  1362.     object *v;
  1363. {
  1364.     return newfloatobject(dgetlongvalue(v));
  1365. }
  1366.  
  1367. static object *
  1368. long_oct(v)
  1369.     object *v;
  1370. {
  1371.     return long_format(v, 8);
  1372. }
  1373.  
  1374. static object *
  1375. long_hex(v)
  1376.     object *v;
  1377. {
  1378.     return long_format(v, 16);
  1379. }
  1380.  
  1381.  
  1382. #define UF (unaryfunc)
  1383. #define BF (binaryfunc)
  1384. #define TF (ternaryfunc)
  1385. #define IF (inquiry)
  1386.  
  1387. static number_methods long_as_number = {
  1388.     BF long_add,    /*nb_add*/
  1389.     BF long_sub,    /*nb_subtract*/
  1390.     BF long_mul,    /*nb_multiply*/
  1391.     BF long_div,    /*nb_divide*/
  1392.     BF long_mod,    /*nb_remainder*/
  1393.     BF long_divmod,    /*nb_divmod*/
  1394.     TF long_pow,    /*nb_power*/
  1395.     UF long_neg,    /*nb_negative*/
  1396.     UF long_pos,    /*tp_positive*/
  1397.     UF long_abs,    /*tp_absolute*/
  1398.     IF long_nonzero,/*tp_nonzero*/
  1399.     UF long_invert,    /*nb_invert*/
  1400.     BF long_lshift,    /*nb_lshift*/
  1401.     BF long_rshift,    /*nb_rshift*/
  1402.     BF long_and,    /*nb_and*/
  1403.     BF long_xor,    /*nb_xor*/
  1404.     BF long_or,    /*nb_or*/
  1405.     (int (*) FPROTO((object **, object **)))
  1406.     (coercion)long_coerce, /*nb_coerce*/
  1407.     UF long_int,    /*nb_int*/
  1408.     UF long_long,    /*nb_long*/
  1409.     UF long_float,    /*nb_float*/
  1410.     UF long_oct,    /*nb_oct*/
  1411.     UF long_hex,    /*nb_hex*/
  1412. };
  1413.  
  1414. typeobject Longtype = {
  1415.     OB_HEAD_INIT(&Typetype)
  1416.     0,
  1417.     "long int",
  1418.     sizeof(longobject) - sizeof(digit),
  1419.     sizeof(digit),
  1420.     (destructor)long_dealloc, /*tp_dealloc*/
  1421.     0,        /*tp_print*/
  1422.     0,        /*tp_getattr*/
  1423.     0,        /*tp_setattr*/
  1424.     (int (*) FPROTO((object *, object *)))
  1425.     (cmpfunc)long_compare, /*tp_compare*/
  1426.     (reprfunc)long_repr, /*tp_repr*/
  1427.     &long_as_number,/*tp_as_number*/
  1428.     0,        /*tp_as_sequence*/
  1429.     0,        /*tp_as_mapping*/
  1430.     (long (*) FPROTO((object *)))
  1431.     (hashfunc)long_hash, /*tp_hash*/
  1432. };
  1433.